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Abstract: Explaining the well established observation that the expansion rate of the uni- 
verse is apparently accelerating is one of the defining scientific problems of our age. Within 
the standard model of cosmology, the repulsive 'dark energy' supposedly responsible has 
no explanation at a fundamental level, despite many varied attempts. A further important 
dilemma in the standard model is the lithium problem, which is the substantial mismatch 
between the theoretical prediction for 7 Li from Big Bang Nucleosynthesis and the value 
that we observe today. This observation is one of the very few we have from along our 
past worldline as opposed to our past lightcone. By releasing the untested assumption 
that the universe is homogeneous on very large scales, both apparent acceleration and the 
lithium problem can be easily accounted for as different aspects of cosmic inhomogeneity, 
without causing problems for other cosmological phenomena such as the cosmic microwave 
background. We illustrate this in the context of a void model. 
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1. Introduction: 

Distant supernovae la (SNIa) explosions are too dim to be explained when taken with other 
cosmological observations, a fact based on a model where the curvature of the universe is 
constant in space and all material in the universe is gravitationally attractive. The standard 
interpretation of SNIa observations within the concordance model of cosmology is that the 
cosmological constant, A, takes on about the same value as the total present day energy 
densities of cold dark matter (CDM) and baryons in the Universe. This extraordinary 
coincidence of numbers, at odds with quantum field theory estimates for the vacuum energy 
by 120 orders of magnitude, leaves many cosmologists slightly queasy. Attempts to explain 
such values typically postulate the physical existence of an infinite number of universes 
with randomly chosen constants of nature, whereby it just so happens that we live in a 
universe like ours at the particular time when A starts to dominate. An elegant idea, but 
untestable except in highly specialised configurations. Alternative explanations for the 
observed brightnesses of SNIa within the homogeneous Friedmann-Lemaitre-Robertson- 
Walker (FLRW) models ether postulate exotic forms of matter with huge negative pressures 
at low temperature, or make ad hoc changes to Einstein's General Relativity on Hubble 
scales [1]. 

A seemingly radical explanation for the SNIa observations suggests that the universe 
is not exactly homogeneous on Hubble scales, but that instead there is significant spa- 
tial variation in the matter distribution, which is accompanied by spatial changes in the 
curvature and local expansion rate. In the simplest of these models, we live in a region 
where the matter density is significantly less than the density of that of the universe on 
super-Hubble scales [2]. The simplicity of the models means that they require violation of 
the Copernican Principle, that we are in some sense 'typical' observers. Consequently, they 
suffer a spatial coincidence problem rather than the temporal one of concordance cosmol- 
ogy. Nevertheless, these models fit SNIa and other local background observations [3, 4, 5], 
and, as we show, the cosmic microwave background (CMB). Structure formation remains 
the unexplored area where these models may fail standard tests, but this is a technically 
challenging problem which has not been properly attempted. The early universe is a key 
area to be analysed, but this is often assumed trivial because the models can evolve from 
FLRW. It's not. 1 

2. Big-Bang nucleosynthesis 

Big-Bang nucleosynthesis (BBN) is the most robust probe of the first instants of the post- 
inflationary Universe. After three minutes, the lightest nuclei (mainly D, 3 He, 4 He, and 
7 Li) were synthesised in observationally significant abundances [7, 8]. Observations of 
these abundances provide powerful constraints on the primordial baryon-to-photon ratio 

1 These models are often mistakenly dismissed because of a so-called weak singularity the models exhibit 
when made to reproduce a negative deceleration parameter at the origin [6]. It is not relevant for cosmo- 
logical modelling because we don't measure a negative deceleration parameter directly, we infer it from a 
parameterised FLRW model. It doesn't occur if we fit a void model to data directly. 
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r] = rib/riry, which is constant in time during adiabatic expansion. In the ACDM model, the 
CMB constrains f] C MB = 6.226 ± 0.17 x 1CT 10 [8] at a redshift z ~ 1100. Observations of 
high redshift low metallicity quasar absorbers tells us D/H= (2.8 ±0.2) x 10 -5 [9] at z ~ 3, 
which in standard BBN leads to r\D = (5.8±0.3) x 10 -10 , in good agreement with the CMB 
constraint. In contrast to these distant measurements at z ~ 10 3 and z ~ 3, primordial 
abundances at z = are either very uncertain (D and 3 He), not a very sensitive baryometer 
( 4 He), or, most importantly, in significant disagreement with these measurements - 7 Li. 
To probe the BBN yield of 7 Li, observations have concentrated on old metal-poor stars in 
the Galactic halo or in Galactic globular clusters. The ratio between rju derived from 7 Li 
at z = and ij£> derived from D at z ~ 3 is found to be r]£)/r]Li ~ 1-5. Within the standard 
model of cosmology, this anomalously low value for i]n disagrees with the CMB derived 
value by up to 5-cr [10]. 

Several attempts at explaining this discrepancy have been done relying on systematic 
errors in the derived abundance (mostly related to parameters of the stellar atmospheres 
such as, e.g., the temperature scale), and uncertainty in the estimates of nuclear reaction 
rates. On the other hand, such corrections are typically found to be much smaller than 
the gap between predictions in standard cosmology and observations (see, e.g., Refs. [11] 
and [10]). Another hypothesis is that 7 Li has been depleted even in metal-poor stars. A 
strong initial astration (i.e., "uniform" depletion affecting the whole Milky Way or even the 
whole Universe) appears rather unlikely (e.g., [12]). Depletion in the stellar atmospheres 
seems instead more feasible. In the globular cluster NGC 6397 [13] and in a sample of 
88 metal poor stars [14], it has been found that describing the depletion through atomic 
diffusion processes partially compensated by turbulence [15] can lead to a Li abundance 
in concordance with cosmological predictions. However, further confirmations are required 
since this mechanism largely depends on a few fine tuned free parameters (see, e.g., the 
discussion in [16] and [17]), depletion by diffusion/turbulence in globular clusters has been 
questioned [18], and the discrepancy still remains in field metal-poor dwarfs when the same 
"depletion correction" of Ref. [13] is applied [12]. Therefore, the state of the art in the 
field is that there is no widely accepted robust astrophysical solution to the 7 Li problem. 

On the other hand, a local value of rj ~ 4 — 5 x 10~ 10 is consistent with all the 
measurements of primordial abundances at z = (see top left panel in Fig. 1). The 
disagreement with high-redshift CMB and D data (probing i] at large distances) shows 
up only when rj is assumed to be homogeneous on super-Hubble scales, as in standard 
cosmology. An inhomogeneous radial profile for ?] can thus solve the 7 Li problem, as we 
show in Fig. 1. 

The local primordial abundances of 7 Li [10], 4 He [19] 2 and 3 He [7] agree with each other 
if rjiQ ~ 4.5 (which we consider as our benchmark value for local rj in the rest of the paper). 

2 For the 4 He mass fraction, we consider Y = 0.2561 ± 0.0108, recently derived in [19], as reference 
value. Previous measurements pointed towards lower values, e.g., Y = 0.240 ± 0.006 reported in [7] which 
is compatible as well with our benchmark scenario. On the other hand, another recent analysis [20] found 

Y — 0.2565 ±0.001 ±0.005, namely, a central value similar to [19], but with significantly smaller error bars. 
This would be in tension with our model as well as with the standard cosmological model, and we consider 

Y from [19] as a more conservative estimate. 
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Figure 1: Constraints on rj. Top left we estimate current constraints on 7710 = 10 w r] from 
different observations. Constraints from 7 Li observations [10] in Galactic globular clusters and 
Galactic halo are shown separately, alongside 4 He [19] and 3 He [7]. These agree with each other if 
7710 ~ 4.5. Local measurements of D show a significant scatter [8] though they also require rjio ~ 4.5 
if we consider the value from Baycsian analyses in [21] and an astration factor of / ~ 2. On the 
other hand, D observations at high redshift (red) [9] and CMB require 7710 — 6. Bottom left we 
show how a varying radial profile for rjio (from ~ 4.5 at the center to ~ 6 asymptotically) can fit 
all the observational constraints, for differing inhomogeneity scales. (The D and CMB constraints 
are in redshift, so move when given in terms of comoving distance r, since r(z) is dependent on 
inhomogeneity profile.) On the right we show the nuclei abundances as a function of z in our 
example model (see text). Filling in points on this graph will test this theory. 



In this scenario, local measurements of D in the interstellar medium (ISM) [21] require in 
turn an astration factor (i.e. depletion of D due to star formation in our Galaxy) of / ~ 2. 
This value lies slightly above the range 1.4 < / < 1.8 predicted in the simulation of [22]. 
However, given the number of assumptions involved in a theoretical model of Galactic 
chemical evolution (see, e.g., the discussion in [23, 24] and reference therein), there is no 
strong argument to exclude / ~ 2 (conservatively, the range / = 1 — 3 can satisfy current 
constraints; as an example see [25]). Moreover, observations of D abundances in the ISM 
show a significant scatter [8]; for sake of concreteness, we reported the value from the 
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analysis in Ref. [21], from which the astration factor follows. Considering a slightly higher 
estimate for the mean local value of D/H leads to smaller / (cf, Fig. 1), while remaining 
consistent with 'Li observations. 

A Gpc-scale inhomogeneity in the baryon-to-photon ratio can therefore explain the 
mismatch between r/n and ijcmb, without violating other constraints on nuclei abundances. 
It should be possible to construct a model of such a variation of 77 within the standard model 
using an O(l) isocurvature mode at early times. However, the fact that such a mode must 
be effectively localised about us (if it wasn't homogeneous at CMB distances it would 
leave a significant distortion in the CMB), and would generate curvature perturbations at 
late times when it enters the Hubble radius (curvature and isocurvature modes are only 
independent on super-Hubble scales), suggests a deeper connection with models that can 
explain supernovae distances without dark energy. 



3. Void models with inhomogeneous 77 

How could dark energy be related to the lithium problem? Let us consider this in the 
context of a toy void model which can fit SNIa data, though our discussion is more widely 
applicable. Spherically symmetric void models have a spatially varying curvature, with 
differing values today at the centre and at large distances (where we observe the CMB). 
The scale of the void is of order the Hubble radius today, so is only directly significant 
at late times. When the curvature mode is larger than the horizon size, gradients in 
the curvature across a distance H{a)~ l are negligible; significant effects appear between 
widely separated radii. For the BBN and local physics at decoupling in CMB calculations 
we assume that early time evolution in a causal patch follows the usual FLRW evolution, 
with a different FLRW model at each radius. 

At early times during the radiation era, we find that unless 77 is fine-tuned to be exactly 
spatially constant, the radiation inhomogeneity grows while the matter inhomogeneity is 
frozen at a value S rn ~ 1 — r/™-' /?^ out ) ^> S r (where (in) refers to quantities at the center 
and (out) at asymptotic distances, and Si = 1 — p$ / p[° with i =matter or radiation, 
and we assume for simplicity that the baryon fraction /& is spatially constant). This is 
shown in Fig. 2 and is the key to linking the lithium and dark energy problems. 

It is not difficult to derive explicitly the early time evolution of the inhomogeneity. At 
early times curvature is negligible and p m <C p T . One can evaluate Si at fixed cosmic time 
by comparing /f/~( out ) = d In a^ n '/d In a^ out ' which, neglecting spatial gradients, leads 
to 
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to first-order in aQ, m /Q, r = p m /Pr, where the fi's always represent todays values. Then, 
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Figure 2: The evolution of inhomogeneity, for an asymptotically flat version of the example 
model. If we evolve an inhomogencous model back to BBN time any inhomogeneity 6 r = 1 — 
p( m ) I p(° ut ) m the radiation component today decays (top left), so that the radiation is homogeneous 
at early times during the radiation era. Meanwhile, the inhomogeneity 5 m present in the matter 
content decays as we head back in time to matter-radiation equality and then grows and freezes. 
This freezing accounts for the 7 Li abundances. On the upper right we show the effect of changing Sr/: 
if this is small the early time matter inhomogeneity S m follows the radiation. On the bottom left, 
we show the evolution of fif and f^ "*-*. The zoomed plot shows the radiation profile changing 
from an overdensity to an underdensity at late times. The evolutions of p( ln ' and p( out ) follow 
from standard FLRW evolutions since spatial gradients are negligible at the center and at large 
distances. Corrections to this assumption would only show up in the dashed line in the blow-up, 
after curvature matter equality. 



integrating and dropping a decaying mode, 
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and using the fact that along a particular worldline 
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we finally have 
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From this we see the results we present graphically, which are presented for the asymptoti- 
cally flat case with /{, = constant: the radiation inhomogeneity is growing linearly with the 
scale factor, while the matter inhomogeneity is frozen very early, decaying slightly from an 
initial value (1 — ) . In the limit r/ m ) = r/ out ) , assumed in the standard model, all this 



rj( out ) 

behaviour disappears, and inhomogeneity grows at the next order in this series expansion; 
in particular, the matter and radiation inhomogeneities grow at the same rate. 

Let us discuss Fig. 2 in more detail. Moving forwards in time, start with an under- 
density in the central matter compensated by a tiny over-density in the radiation, required 
for the models considered here. The matter density grows relative to radiation as we 
leave the radiation era, a process which starts to smooth out existing inhomogeneity in 
the matter, while decreasing the rate of growth of the inhomogeneity in the radiation. 
At the centre of the void the curvature starts to become significant at t ~ 1 Gyr at the 
expense of the matter component, but this doesn't happen so fast where the curvature 
starts smaller, at large distances such as at the CMB (where we observe 7/( out )). In the 
language of perturbation theory, the curvature mode enters the Hubble radius and starts 
to grow. The growth of the curvature at the centre amplifies the late time inhomogeneity 
and is the process which accounts for the SNIa magnitudes. The key finding is that the 
difference in the baryon density between the centre of the void and the outside is significant 
and 0(1) at BBN time. While the model homogenises as t — > 0, it does so only in the 
dominant radiation component, and not in the matter. Only when the matter becomes 
relativistic will the model truly homogenise. 

The 0(1) difference in the baryon density inside and outside the void implies an 0(1) 
difference in the spatial profile 77 at BBN - exactly what we measure. Furthermore, SNIa 
observations imply a decrease in the local matter density, the same as Li observations. 
Hence, a void-like profile for rj with our position close to the centre of the void can easily 
give an explanation to the 7 Li problem if ijcmb = ?/ out ) ~ 1.5 x 7y( m ); this model then 
evolves into a void which can explain dark energy. 3 Relevant BBN processes depend only 
on the FLRW parameters for each patch; since we do not consider exotic contributions 



3 A recent analysis [26] constrained large-scale inhomogeneities in 77 using D abundances (observed in 
different directions of the sky) and galaxy cluster gas fraction [27] (observed up to z ~ 1). Because we 
consider a spherically symmetric inhomogeneity and a homogeneous baryon fraction ft = f2&/f2 m (with 
ft, ~ 0.17 from the CMB), these bounds are automatically satisfied. 
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to Hubble rate or neutrino-antineutrino asymmetry, the deviation from standard BBN is 
only given by an inhomogeneous profile for rj giving different abundances in each patch. 4 
The spatial distribution of nuclei primordial abundances today reflects the inhomogeneity 
at BBN time, as there is no homogenisation by structure formation on Gpc-scales. Note 
that we can fine-tune the initial conditions such that the lithium problem is solved, but 
without creating a void at late times (by enforcing a constant-curvature model). This 
requires an 0(1) difference in the radiation density profile today, and a total density profile 
inhomogeneity which grows to ~ 0.05 at equality, then decays. This is much larger than a 
standard density perturbation at equality. Thus, there appears no way to solve the lithium 
problem within the standard FLRW models using this setup, though it might be possible 
within perturbed FLRW models with a large spherical isocurvature mode. 

4. An example 

In the previous Section, we described how models with Hubble-scale inhomogeneity can 
solve the 7 Li problem. Now we show with an example that this is not only a general 
possibility, but rather the actual parameters involved in this scenario are in the same 
ball-park required by other cosmological observations, such as Hubble rate measurements, 
supernovae observations and the CMB. The calculation presented here is discussed in detail 
in [30]. 

Voids may be described by a profile for the present-day matter density f2 m (r) = £% ut ^ — 
(flm Ut ^ — &m ) e~ r2 / 2a2 = 1 — f^fe(r), and a similar profile for the radiation density (usually 
ignored but important here). The Hubble rate H(t,r) is given by a generalized Friedman 
equation (see Appendix in [30]) and we choose a(to, r) = 1 to fix the gauge freedom in void 
models. For our purposes, we can calculate nearly everything along the central worldline 
and along one at the radius of the CMB (r ^> a), neglecting all r-dependence, and use 
the standard Friedman equation. The exception is the area distance to the CMB, which is 
calculated in the normal way (see e.g., [5]), and is affected by the shape and width of the 
void. Indeed CMB fluctuations at small scales (and focusing on the first peaks) are mainly 
set by three parameters [31, 32, 33]: the angular scale of the sound horizon l a , the particle 
horizon at matter-radiation equality l eq , and the baryon-photon density ratio at decoupling 
R* = 3 p~f\^, where l a /l e q and R* depend only on local physics at decoupling, and then 
the area distance to the last scattering surface enters in the same way in l a and l eq (for a 
detailed analysis, see [30]). In the example we are going to present, we enforce l a = 302, 
l e „ = 136, and R* = 0.63, in near perfect agreement with WMAP data [34]. 

We fix the void parameters on the inside and outside using sample observations as 
follows: 

4 Our inhomogeneous BBN is completely different from the inhomogeneous BBN relying on a vary- 
ing baryon to photon and/or neutron to proton ratio on scales comparable to particle diffusion horizon 
at BBN [28]. We also differ from non-standard BBN with large but not giant scale inhomogeneities (or 
anisotropics) [29] (i.e., scales smaller than (or comparable to) astrophysical structure scales), where primor- 
dial abundances produced in different regions are then mixed by structure formation. 
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Local constraints: Take T (in) rj 2.725K => ft^)/^ ~ 2.49 x 10~ 5 from the CMB 

temperature today (which gives 0^ ~ for iV e g = 3.04); r/( m ) ~ 4.5 x 10 -10 from 

7 Li measurements; = 0.74 ± 4 from HST measurements [35]. We expect ' ~ 0.15 
from SNIa constraints. The age is to = Jq ^jj along the central worldline. We assume the 
Bang time to be homogeneous tg = 0, so to is the age of the universe everywhere today. 
Asymptotic constraints: Given CMB observations of l a ,leq and R*, we may fix r]( out ' 
and fff from i?* and l a /l eq , and then find cIa(z*) from l a or l eq . This implies ?/ out ) 6.2x 
10~ 10 , and / 6 ~ 0.17 [30] (in our example, we assume /& =const., so that / 6 = / 6 )• 
We can also require ?/ out ) 6.2 x 10~ 10 to be compatible with D observations at z ~ 3 
(although the void profile and width can affect this, if the asymptotic value is not reached 
by this redshift). The redshift of decoupling is defined by 1 + = T* /Tq \ where 
the local temperature of the CMB, T*, can be in first approximation determined when the 
Thompson scattering interaction rate T = H^ ont \ and making use of the Saha equation, or 
more precisely through a fitting formula [30]. 

Constraints from distances: The matter profile is constrained by measurements of 
distances. The matter density has to be O m = 1 — ~ 0.1 — 0.2 inside the void at 
the centre, which must raise to £l m ~ 0.5 by z ~ 1, or a (co-moving) distance of a few 
Gpc, to satisfy SNIa observations [5, 36] (then may or may not reach asymptotic flatness). 
Our void must then be able to fit the area distance to the last scattering surface (Ia(z*) ~ 
13Mpc derived from CMB observations, which constrains a and fim"''' Consistent pictures 
typically require a ~ few Gpc, and asymptotically open models, although considering 
profiles more complicated than a Gaussian (or a significantly lower h^), > 1 is also 

possible [30]. 

Derived parameters: The baryon density is given by: 0[ in) /i2. n) «3.70xl0 6 r ? ( in )(T (in) /2.7K) 3 . 

This fixes = ^[ m ^ /f^m^ ■ We may estimate /j( out ) from to = j ( Q u^) jj-(out) ' i§ nor i n § 

the contribution from Qr° ut ^ which we don't yet know. This introduces a tiny error, which 
can be accounted for iteratively. The effective FLRW temperature of the CMB today at 
the radius where we see the CMB emitted is: 



^(out) fo {ont) h 2 (in) 



. 1/3 

tp"""' = / " m "(out) 

T (in) \ tffhfa ^ ( ° ut) 



(4.1) 



Note that T (out) is significantly different from T (in) . Then Q<° ut) ^ out) 4.48xl0^ 7 (T (out) /lK) 4 
follows. All 'out' parameters which are set at t = to refer to an asymptotic observer 
(who can also be seen as defining an effective FLRW model) which has the same intrinsic 
CMB (though usually different area distance) as the void observer will see 5 . To calculate 
the full CMB spectrum, calculate the shift S = d L A TB {z = T* (out) /T^ n) )/d^ LRW(ont) (z = 
Tj out ^ /Tq° u ^ ) for the angular power spectrum output from a standard CMB code set to 



5 Mostly, however, they will be close to the real values along the wordline of the CMB (but note that they 
aren't actually required in this analysis). For example, a real outside observer (i.e., along the wordline of the 
CMB) will see a cold spot in the CMB where the void is, which will lower the mean 'outside' temperature. 
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Figure 3: The normalised CMB angular power spectrum. This is for the example model 
described in the text, shown on a log scale for I < 200 and linear for I > 200. The power spectrum is 
shown against a default flat concordance model with zero tilt, and h — 0.7, £l m h 2 = 0.137, f^/i 2 = 
0.0226. There's nothing between the two models for high /, with the maximum difference around 
1%. This is also the case for the polarisation spectra (not shown). 



the 'out' values (importantly using T ) to match the data [37]. This gives CJ md oc 
with any amplitude shift being absorbed by the unknown amplitude in the 
primordial power spectrum. The result is shown in Fig. 3. 

From this procedure one can construct a void model in agreement with CMB, SNIa, 
HST and BBN observations (see Fig. 4) 6 . As an example, we choose h = 0.70, r/ 1 "* 1 = 
4.5 x 10 -10 , fb =const. (which implies ~ 0.2) a Gaussian profile which has a full 

width at half maximum of ~ 3Gpc and has £l m ~ 0.7 at the CMB distance. It provides 
an excellent fit to the CMB, as shown by the angular power spectrum in Fig. 3. It is 
interesting to note that if we make ?/ in ) larger than this, then it makes flffl too high for 
the SNIa (we loose the 'bump' in the distance modulus) - so, in effect the SNIa favour a 
low rj( m \ 7 For a more detailed description of our CMB computation see Ref. [30]. 

This model fits the Li constraints in the most natural way: it requires an 0(1) decrease 
in the baryon density at the centre, just as dark energy requires an 0(1) decrease in the 
CDM. While there still exists freedom in the models to fit more complicated possibilities 
(such as an inhomogeneous bang time or fb) 8 , we do not need to invoke them. (A full 
likelihood analysis will be explored elsewhere.) There is no fine-tuning involved in these 
fits. By contrast, the assumption in the standard model that rj is constant across comoving 
volumes ~ Gpc 3 may be considered extremely fine-tuned in an inhomogeneous context. 

The large scale CMB together with the Baryon Acoustic Oscillations (BAO) depend 

6 Age constraints [38] favour a slightly lower value of h than the HST value, but with Q m low enough to 
satisfy SNIa data, these are relatively easy to satisfy. 

7 Previous analysis of the CMB in void models found that the CMB couldn't be accounted for except in 
very unusual configurations [39, 37], although [37] show that an inhomogeneous bang time allowed a wide 
class of voids to fit. Their analysis is based on test radiation in pure dust void models, and we find that 
including inhomogeneous radiation in the model allows a much wider class of voids to fit the CMB [30]. 

8 Note that they can also be invoked in order to compensate for possible deviations from FLRW evolutions 
we consider at the center and asymptotically. 
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Figure 4: A spacetime overview of an example model. At the three key times we show 
the matter profile (solid). The radiation profile (dashed) at early times is an over density in this 
scenario. 



on the detailed evolution of perturbations during the curvature era which are not yet 
understood in the context of inhomogeneous models. Moreover to extract constraints from 
BAO and kinetic Sunyaev-Zeldovich a detailed treatment of the radiation is necessary [30]. 

While important for CMB-related measurements, the evolution of radiation at late 
times is completely irrelevant for our BBN analysis. For example, suppose one can construct 
models where the temperature of the radiation does not evolve as 1/a along the central 
worldline at late times (i.e., in the curvature era), which leads to non-conservation of 7] 
(evaluating the radiation temperature using the geodesies of the LTB solution leads to this 
conclusion). While in this case the CMB predictions will be different [39, 37], the core 
of our analysis is unchanged because 7 Li abundances constrain at early times, when 
FLRW evolution in each causal patch is certainly a good approximation for a Gpc-scale 
inhomogeneity. In particular, Fig. 2 and surrounding analysis is unchanged except in the 
curvature era. 

5. Conclusions 

We have found that the inhomogeneous models which explain dark energy without exotic 
physics at late times also explain the lithium abundances in a natural way. In retrospect, 
the fact that 7 Li gives us an unusual glimpse of our past worldline implies that this is 
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actually a test of homogeneity at very early times. We have argued that, while there 
may exist an explanation of the lithium abundances in terms of an anomalously large Gpc 
isocurvature mode in FLRW, the fact that this must be spherical or otherwise localised 
to avoid problems with the CMB implies a strong link with the spherical inhomogeneous 
models for dark energy. The origin of the inhomogeneous profile lies in the primordial 
baryon asymmetry- entropy density ratio, and may in fact be laid down at the onset of 
an inflationary period which lasts slightly less than the 60-efolds or so required for the 
standard model; it cannot plausibly exist within standard slow roll inflation. At first sight 
the void models appear fine-tuned as we must be within tens of Mpc of the centre [40, 41], 
a coincidence of 1 part in ~ 10 8 . But this has to be compared to coincidences of 10 120 
or even 10 500 encountered in the standard model. In any case, the spherically symmetric 
void models are merely a simple first attempt at an inhomogeneous universe and so should 
be understood as toy models in this context. We would argue that more sophisticated 
models of inhomogeneity should be explored to address the Copernican problem, and that 
the problems of structure formation and inflation should now be properly investigated. 

Acknowledgments 

We thank Bruce Bassett, Timothy Clifton, George Ellis, Pedro Ferreira, Roy Maartens and 
Jean-Philippe Uzan for discussions. This work is supported by the National Research Foun- 
dation (South Africa). M.R. acknowledges funding by the Centre for High Performance 
Computing, Cape Town and the hospitality of the African Institute for Mathematical Sci- 
ences. 

Note added 

Since the first version of this paper appeared on the arXiv, a number of other papers 
have appeared which also discuss the CMB in void models [42, 43, 44, 45] (see also [46]). 
Comparisons with our work are discussed in our follow-up paper [30]. 
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